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INTRODUCTION 

Colossal magnetoresistance in metallic oxides such as Ri_ x X x Mn03 (where R = La, 
Pr, Nd; X = Sr, Ca, Ba, Pb) is attracting considerable attention |j] due to its potential 
technological applications. A variety of experiments have revealed that oxide mangan- 
ites have a rich phase diagram [0 with regions of antiferromagnetic (AF) and ferromag- 
netic (FM) order, as well as charge ordering, and a peculiar insulating state above the 
FM critical temperature T c . Recently, layered manganite compounds Lai.2Sri. 8 Mn 2 07 
have also been synthesized || with properties similar to those of their 3D counterparts. 
Strong correlations are important for transition-metal oxides, and, thus, theoretical 
guidance is needed to understand the behavior of manganites and for the design of new 
experiments. 

The appearance of ferromagnetism at low temperatures can be explained using the 
so-called Double Exchange (DE) mechanism j|, ||. However, the DE model is incom- 
plete to describe the entire phase diagram observed experimentally. For instance, the 
electron-phonon coupling may be crucial to account for the insulating properties above 
T c |J. The presence of a Berry phase at large Hund-coupling also challenges predic- 
tions from the DE model 0. In a series of recent papers (Refs. || §, |K], PH ) we have 
remarked that another phenomena occurring in manganites and not included in the 
DE description, namely the charge ordering effect, may be contained in a more funda- 
mental Kondo-like model. More precisely, in those papers ^ |], |TD[ it was reported the 
presence of phase separation (PS) between hole-undoped antiferromagnetic and hole- 
rich ferromagnetic regions in the low temperature phase diagram of the one-orbital FM 
Kondo model. In a recent paper by the same authors a similar phenomenon was also 



observed using two orbitals and Jahn- Teller phonons [11 . In this case the two phases 



involved are both spin ferromagnetic and the orbital degrees of freedom are responsible 
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for the phase separation. Upon the inclusion of long-range Coulombic repulsion, charge 
ordering in the form of non-trivial extended structures (such as stripes) could be sta- 



bilized, similarly as discussed for the cuprates [12, n3L 041 O but now also including 



ferromagnetic domains. The analysis carried out by our group suggests that phenom- 
ena as rich as observed in the high-Tc superconductors may exist in the manganites as 
well, and hopefully our effort will induce further theoretical and experimental work in 
this context. 

The present paper should be considered as an informal review of the present status 
of the computational studies of models for manganites that reported the existence of 
phase separation at low temperatures. It does not pretend to be a comprehensive 
article, and thus we encourage the readers to consult the literature mentioned here to 
find out additional references. The paper is organized as follows: in the next section 
the results for the case of the one-orbital Kondo model are reviewed, with emphasis 
on the phase diagram. Most of the calculations are performed with classical t 2g -spins 
but results for quantum spins are also shown. In addition, models that interpolate 
between Cu-oxides and Mn-oxides have also been studied, as well as the influence on 
our conclusions of Coulomb interactions beyond the on-site term. In the following 
section results recently reported for the case of the two-orbital model with Jahn- Teller 
phonons are described. Once again the main emphasis is given to the phase diagram. 
In both of those previous sections the main result is that tendencies to ferromagnetism 
and phase-separation are in strong competition in these models. Such a phenomenon 
appears so clearly in all dimensions of interest and for such a wide range of models 
that it leads us to believe that it may be relevant for experiments in the context of 
manganites. A discussion of experimental literature that have reported some sort of 
charge inhomogeneity compatible with phase separation (after Coulomb interactions 
are added) are briefly described towards the end of the paper. In the last section a 
summary is provided. 



ONE ORBITAL FERROMAGNETIC KONDO MODEL 
Results for Classical t 2g -Spins 

The FM Kondo Hamiltonian M, is defined as 



H = -tJ2 (4t c J<t + h.c.) -JhY, c L^Ci/? • Si, (1) 

(ij)o- ia/3 

where c\ a are destruction operators for one species of e 9 -fermions at site i with spin a, 
and Si is the total spin of the t 2g electrons, assumed localized. The first term is the e g 
electron transfer between nearest-neighbor Mn-ions, Jh > is the Hund coupling, the 
number of sites is L, and the rest of the notation is standard. The density is adjusted 
using a chemical potential //. In this section and Refs. || |9j the spin Si is considered 
classical (with |Si| = 1), unless otherwise stated. Although models beyond Eq.(l) 
may be needed to fully understand the manganites (notably those that include lattice 
effects as studied later in this paper), it is important to analyze in detail the properties 
of simple Hamiltonians to clarify if part of the experimental rich phase diagram can be 
accounted for using purely electronic models. 

To study Eq.(l) in the ^-spin classical limit a Monte Carlo (MC) technique was 
used: first, the trace over the e 9 -fermions in the partition function was carried out 
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exactly diagonalizing the 2L x 2L Hamiltonian of electrons in the background of the 
spins {Si}, using library subroutines. The fermionic trace is a positive function of the 
classical spins and the resulting integration over the two angles per site parametrizing 
the Si variables can be performed with a standard MC algorithm without "sign prob- 
lems" . In addition, part of the calculations of Refs. || |J were also performed with the 
Dynamical Mean-Field approximation (D=oo) [[L(|, the Density-Matrix Renormaliza- 



tion Group (DMRG), and the Lanczos method. Special care must be taken with the 
boundary conditions (BC) |T7], IE . 

Our results are summarized in the phase diagram of Fig.l. In ID (and also in 2D, 
not shown) and at low temperatures clear indications of (i) ferromagnetism, (ii) spin 
incommensurate (IC) correlations, and (iii) phase separation were identified. Results 
are also available in small 3D clusters and qualitatively they agree with those in Fig.l. 
The same occurs working at D = oo (see Refs. || fj). In ID we also obtained results 
with quantum t 29 -spins S=3/2 (shown below) which are in good agreement with Fig.l. 




Figure 1: Phase diagram of the FM Kondo model reported in Refs. [8, 9] obtained 
using numerical techniques. FM, IC, and PS denote regimes with FM correlations, 
spin incommensurate correlations, and with phase separation between undoped AF 
and hole-rich FM regions, respectively. The results shown correspond to one dimension 
and they were obtained with MC simulations at T = t/75 using chains with L=20, 30 
and 40 sites. This result is taken from Refs. [8, 9]. 

The boundaries of the FM region of the phase diagram were found evaluating the 
spin-spin correlation defined as S(q) = (1/L) Z)j m e l ^~ m "> ' q (Sj • S m ). Fig.2 shows S(q) 
at zero momentum vs. T/t for typical examples in ID and 2D. The rapid increase of the 
spin correlations as T is reduced and as the lattice size grows clearly points towards the 



existence of ferromagnetism in the system [jl9|j . It is natural to assume that the driving 



force for FM in this context is the DE mechanism. Repeating this procedure for a variety 
of couplings and densities, the robust region of FM shown in Fig.l was determined. In 
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the small Jn/t region IC correlations were observed monitoring S(q) Both in ID 
and 2D there is one dominant peak which moves away from the AF location at (n) = 1 
towards zero momentum as (n) decreases. In the 2D clusters the peak moves from 
(ir, 7r) towards (vr, 0) and (0, ir), rather than along the main diagonal. Note that our 
computational study predicts IC correlations only in the small and intermediate Jn/t 
regime. 
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Figure 2: Spin-spin correlations of the classical spins at zero momentum 5(q = 0) 
vs temperature T (in units of t). MC results for several lattice sizes are shown on (a) 
chains and (b) 2D clusters. The density and coupling are shown. In (a) closed shell 
BC are used i.e. periodic BC for L = 10 and 30 and antiperiodic BC for L = 20 and 
40. In (b) open BC are used. Results taken from Ref. [8]. 

The main result of Refs. [§, |9[] is contained in Fig. 3 where the computational evidence 
for the existence of phase separation in dimensions 1,2, and oo is given. The presence of 
a discontinuity in (n) vs fj, shows that some electronic densities can not be stabilized by 



tuning the chemical potential |2T] . If the system is nominally prepared with such density 
it will spontaneously separate into two regions with the densities corresponding to the 
extremes of the discontinuities of Fig.3. By analyzing these extremes the properties 
of the two domains can be studied. One region is undoped ((n) = 1) with strong 
AF correlations, while the other contains all the holes and the spin-spin correlations 
between the classical spins are FM (see the inset of Fig. 3a. The results are similar in 
D=2 and infinite). This is natural since holes optimize their kinetic energy in a FM 
background. On the other hand, at (n) = 1 the DE mechanism is not operative: if the 
electrons fully align their spins they cannot move in the conduction band due to the 
Pauli principle. Then, energetically an AF pattern is formed. As Jh grows, the jump 
in Fig.3 is reduced and it tends to disappear in the Jh = oo limit. 

Experimentally, PS may be detectable using neutron diffraction techniques, since 
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S(q) should present a two peak structure, one located at the AF position and the other 
at zero momentum. Since this also occurs in a canted ferromagnetic state care must be 
taken in the analysis of the experimental data. Another alternative is that Coulombic 
forces prevent the macroscopic accumulation of charge intrinsic of a PS regime. Thus, 
microscopic hole-rich domains may develop in the manganites, perhaps in the form of 
stripes. 
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Figure 3: Electronic density (n) vs chemical potential fi in (a) D=l, (b) D=2, and 
(c) D = oo clusters. Temperatures are indicated. The coupling is Jnjt = 8.0 in (a) 
and (b) and Jh/W = 4.0 in (c) (for the definition of W see Ref. [8]). PBC were used 
both in D=l and 2. The discontinuities shown in the figures are indicative of PS. In 
(a) the inset contains the spin-spin correlation in real space at densities 1.0 and 0.72 
showing that indeed PS occurs between AF and FM regions. The results are taken 
from Refs. [8, 9] where more details can be found. 

Although the phase diagram of Fig.l has PS close to half-filling, actually this phe- 
nomenon also occurs at (n) ~ if an extra direct AF exchange interaction between the 
localized spins is included |]10||. This coupling may be originated in a small hopping 



amplitude for the t 2g electrons. At (n) = 0, model Eq.(l) supplemented by a Heisen- 
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berg coupling J' ft among the localized spins produces an AF phase, as in experiments, 
which upon electron doping induces a competition between AF (with no e g -electrons) 
and FM electron-rich regions, similarly as for (n) = 1 but replacing holes by electrons. 
Thus, PS or charge ordering could exist in manganites both at large and small fermionic 
densities. A careful study of the influence of the J' coupling on the phase diagram of 



the 1-orbital Kondo model has been recently presented in Ref. JlO 

For completeness, upper bounds on the 3D critical temperature T c were also pro- 
vided in Refs. || [|. Using MC simulations in principle it is possible to calculate T c 
accurately. However, the algorithm used in our studies prevented us from studying 
clusters larger than 6 3 even at Jh — oo. In spite of this limitation, monitoring the 
spin-spin correlations in real space allows us to judge at what temperature T* the cor- 
relation length reaches the boundary of the 6 3 cluster. Since the bulk T c is smaller than 
T*, this give us upper bounds for the critical temperature. Following this procedure at 
(n) = 0.5, it was found that for T~0.lt robust correlations reach the boundary, while 
for T > 0.12t the correlation is short-ranged. Thus, at this density we estimate that 
T c < 0.12t. Since results for the e g electrons bandwidth range from BW ~ 1 eV |22 



to BW ~ 4eV |23|], producing a hopping t = BW/ 12 between 0.08 and 0.33 eV, 



then the estimations for the critical temperature of Refs. f||, |J range roughly between 
T c ~ 100 K and 400 K. This is within the experimental range and in agreement with 



other results [16|, 24j| . Then, in principle purely electronic models can account for T c . 



Quantum t 2g -Spins: Phase Diagram 

In the previous section we have used classical localized spins to simplify the technical 
aspects of the computational study and since the realistic values of those spins is large 
(3/2). However, it is interesting to study the phase diagram for the case of quantum 
t 2g -spins at least in one special case, in order to verify that indeed the use of classical 
spins is a good approximation. Using the Lanczos and DMRG techniques, the phase 
diagram for the cases of t 2g -spins S = 3/2 and 1/2 was reported in Refs. || || for 
the case of one dimension, and one of them is shown in Fig. 4. The techniques used 
work in the canonical ensemble, instead of grand-canonical as the MC simulation of 
the previous section. The presence of phase separation is investigated here using the 
compressibility (which if it becomes negative signals the instability of the density under 
study). Fig. 4 shows that the three regimes found before (namely PS, FM, and IC) are 
also observed in the quantum case. The results for 5 = 3/2 are even quantitatively 
similar to those found with the MC technique of the previous section, while those for 
5=1/2 are only qualitatively similar. The overall conclusion is that the use of classical 
t 2g -spins in the Ferromagnetic Kondo Model was found to be a good approximation, 
and no important differences in the phase diagram have been observed (at least in this 
one-dimensional example) when quantum spins are used. Studies in dimensions larger 
than one would be too involved using quantum localized spins, and thus the similarities 
between the classical and quantum cases found here help us to investigate those more 
realistic dimensions. 



Models for Other Transition-Metal Oxides 



In the previous sections we have found that models for manganites analyzed with 
reliable unbiased techniques provide a phase diagram where there are two main phases 
in competition: a ferromagnetic regime and a "phase-separated" regime (actually the 
latter is not a phase but a region of unstable densities). Such a clear result suggest that 
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Figure 4: Phase diagram of the Ferromagnetic Kondo model using quantum t 2g -spins. 
The results correspond to a localized spin of value S = 3/2 (normalized as \S\ = 1, 
see Ref. [9]). The notation is the same as in Fig.l. The techniques used to obtain the 
results are Lanczos and DMRG. The result is taken from Ref. [9]. 



the phenomenon maybe more general than expected. Actually recent studies of models 



for transition- metal oxides in general [|17]], in the limit of a large Hund-coupling, revealed 
phase diagrams where indeed FM and PS dominated. Let us consider the evolution 
of such a phase diagram as we move in the periodic table of elements from copper, to 



nickel, etc, eventually arriving to manganese. The study of Ref. |17j was carried out in 
one-dimension for technical reasons, but nevertheless the tendencies observed are strong 
enough to be valid in other dimensions as well. For the case of copper the model used 
was the plain "t-J" model, namely a model where we have spins of value 5=1/2 and 
holes which carry no spin. The phase diagram in this case was known for a long time. 
It contains no indications of ferromagnetism (at least no "clear" indications), and there 



is phase separation only at large J/t. However, the study of Ref. |T7j for Ni-oxides using 
a t-J-like model with "spins" of value S — 1 and "holes" with spin S = 1/2 (since they 
correspond to the removal of one electron of a Ni-site) revealed that both FM and PS 
are clearly enhanced in the phase diagram with respect to the case of Cu-oxides. The 
technical details can be found in Ref. [17]. As the spin grows these tendencies become 



stronger and actually for the case of the Mn-oxides only a small window separates the 
fully polarized FM and PS regimes. This result suggests that the phenomenon is very 
general and should be searched for in other transition-metal oxides as well. 

Influence of Coulomb Interactions 

The accumulation of charge that occurs in the phase separation regime described 
in the previous section is not stable once Coulombic interactions beyond the on-site 
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term are added to the problem. It is expected that in the presence of these long- 
range interactions the ground state will have clusters of one phase embedded into 
the other (droplet regime). However, a computational study including these Coulomb 
interactions is very difficult. For instance, the analysis of model Eq.(l) would require 
a Hubbard- Stratonovich decomposition of the fermionic interaction which not only 
introduces more degrees of freedom, but in addition the infamous "sign problem" which 
will prevent the study from being carried out at low enough temperatures. Then, 
only techniques such as Lanczos or DMRG can handle properly this type of models. 
Recently, the analysis including Coulombic terms started by using the "minimum" 
possible Hamiltonian, namely the 1-orbital Ferromagnetic Kondo Model with localized 
spin-1/2 and with the addition of an on-site Coulomb term and a nearest-neighbor 
(NN) density-density l^-coupling that penalizes the accumulation of charge in two NN 
sites 1 . The results are still being analyzed so here only a qualitative summary will 



be presented. 

The conclusions of the study of Malvezzi et al ||25|| . thus far are the following: (1) 



the on-site coupling U apparently does not affect qualitatively the results found at 
U = as long as the Hund-coupling is large; (2) in the regime of phase separation 
near (n) = 1 the addition of a V^-term produces results compatible with the cluster 
formation anticipated in the previous paragraph. Actually even at finite V in a regime 
where the compressibility shows that phase separation no longer exists (i.e. all densities 
are stable), clear remnants of such a PS regime are found. For instance, at V — 
the spin structure factor S(q) has robust peaks at both q = and tt. Introducing 
V, these peaks still exist but they move slightly from the and tc positions forming 
incommensurate structures; (3) in the regime where at V = there is ferromagnetism 
induced by a double-exchange process, the inclusion of a nonzero coupling V transforms 
the ground state into a FM charge- density- wave (CDW). In this phase the spins are 
aligned but the charge is not uniformly distributed. The overall phase diagram is rather 
complicated. We are actively working in this context to analyze the qualitative aspects 
of such a result. But it is already clear that charge-ordering tendencies are observed at 
all densities upon introducing V, in agreement with experiments. More detailed results 
will be presented soon. 



TWO ORBITALS AND JAHN-TELLER PHONONS 
Phase Diagram 

In spite of the rich phase diagram observed in the study of the 1-orbital Ferromag- 
netic Kondo Model in the previous section and its similarities with experiments, there 
are aspects of the real manganites that require a more sophisticated approach. For 
instance, dynamical Jahn- Teller (JT) distortions are claimed to be very important 0, 
and a proper description of the recently observed orbital order |26j obviously needs at 
least two orbitals. Such a multi-orbital model with JT phonons is nontrivial, and thus 
far it has been studied only using the dynamical mean- field approximation ||27|| . The 



previous experience with the 1-orbital case suggests that a computational analysis is 
actually crucial to understand its properties. In addition, it is conceptually interesting 
to analyze whether the PS described before || [TO] exists also in multi-orbital models. 

The first computational study of a 2-orbital model for manganites including JT 
phonons was reported recently by the authors in Ref. [[□]] and here a summary of the 
main conclusions will be discussed. The results show a rich phase diagram including 



S 



a novel regime of PS induced by the orbital, rather than the spin, degrees of freedom 
(DOF). The Hamiltonian used in that study had three contributions H^jt — Hk + 
Hjt + Haf- The first term is 

H K = - ^( c L c J6<t + h.c.) -J H J2Si- c\ aa a a pc iaP , (2) 

(ij)crab iaa/3 

where (ij) denotes nearest-neighbor lattice sites, Jh > is the Hund coupling, the 
hopping amplitudes t ab are described in Ref. p8| , a, b = 1,2 are the two e 9 -orbitals, 
the t2 9 spins Si are assumed to be classical (with |Si| = 1) since their actual value in 
Mn-oxides (3/2) is large |p9| |, and the rest of the notation is standard. None of the 



results described below depends crucially on the set {t a b} selected P8| . The energy 
units are chosen such that t\\ = 1 in the x-direction. In addition, since Jh is large 
in the real manganites, here it will be fixed to 8 (largest scale in the problem) unless 
otherwise stated. The e 9 -density (n) is adjusted using a chemical potential /i. 
The coupling with JT-phonons is through || [3(| 

Hjt = A £ 4 aa Q?c iba + \ £(Q (2) ? + Q^), (3) 

iaba i 

where Q\ l = —Qf 2 = Q[ , and Q\ 2 = Q 21 = qI . These phonons are assumed to 
be classical, which substantially simplifies the computational study. This is a reason- 
able first approximation towards the determination of the phase diagram of the H^jt 
model. Finally, a small coupling between the ^-spins is needed to account for the 
AF character of the real materials even when all La is replaced by Ca or Sr (fully 
doped manganites). This classical Heisenberg term is Haf = J'J2{ij) Si • Sj, where J' 



is fixed to 0.05 throughout the paper, a value compatible with experiments f3~Ifl . To 
study H KJT a Monte Carlo (MC) technique similar to that employed in Refs. |S|, [| [LTJ 
and in the previous section for the 1-orbital problem was used. Finally, to analyze 
orbital correlations the pseudopin operator Ti = | J^aab Aao° ahC-Voo was used, while for 
spin correlations the operator is standard. The Fourier-transform of the pseudospin 
correlations is defined as T(q) = j Xa,m e * q (T m ■ Ti), with a similar definition for 
the spin structure factor S(q). 

Let us first consider the limit (n) = 1.0, corresponding to undoped manganites. 
Fig. 5 shows T(q) and S(q) at representative momenta q = and 7r vs. A. For small 
electron-phonon coupling the results are similar to those at A = 0.0, namely a large 
S(0) indicates a tendency to spin-FM order induced by DE (as in the qualitatively 
similar 1-orbital problem at (n) = 0.5 ||). The small values of T(q) imply that in 
this regime the orbitals remain disordered. When the coupling reaches A c i ~ 1.0, the 
rapid increase of T(ir) now suggests that the ground state has a tendency to form a 
staggered (or "antiferro" AF) orbital pattern, with the spins remaining FM aligned 
since 5(0) is large. The existence of this phase was discussed before, but using multi- 



orbital Hubbard models with Coulomb interactions and without phonons p2| . Our 
results show that it can be induced by JT phonons as well. As the coupling increases 
further, another transition at A C 2 ~ 2.0 occurs to a spin-AF orbital-FM state {S(n) and 
T(0) are large). In this region a 1-orbital approximation is suitable. 

The three regimes of Fig. 5 can be understood in the limit where A and Jh are the 
largest scales, and using t 12 = t 2 \ = 0, t u = t 2 2 = t for simplicity. For parallel spins 
with orbitals split in a staggered (uniform) pattern, the energy per site at lowest order 
in t is ~ —t 2 /A (~ 0), where A is the orbital splitting. For antiparallel spins with 
uniform (staggered) orbital splitting, the energy is ~ —t 2 /2Jn (~ —t 2 /(2Jn + A)). 
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Figure 5: (a) T(q) and S(q) vs A, working at (n) = 1.0, T = 1/75, Jh = 8, J' = 0.05, 
and in ID with 10 sites. {t a fc} correspond to set Ti (see [28]). Results with sets T 2 
and T 3 are qualitatively the same; (b) Same as (a) but for a 4 2 cluster, T = 1/50, and 
hopping T 3 (T 4 ) in the y (x) direction, q = 0(h) denotes (0,0) ((ir, 7r)); (c) Same as (a) 
but for a 4 3 cluster, T = 1/50, the 3D hopping amplitudes of Ref. [28], and Jh = oo. 
q = 0(7r) denotes (0,0,0) ((7r, tt, tt)). These results were taken from Ref. [11]. 



Then, when A < 2Jh ("intermediate" As), a spin-FM orbital- AF order dominates, 
while as A grows further a transition to a spin-AF orbital-FM ground state is expected. 
Note that this reasoning is actually valid in any dimension. To confirm this prediction, 
results for a 2D cluster in Fig. 5b are shown. Indeed the qualitative behavior is very 
similar to the ID results. In 3D (Fig.5c), where studies on 4 3 clusters can only be done 
at Jh = oo to reduce the number of DOF, at least two of the regimes of Fig. 5a have 



been identified 33 



Let us discuss the transport properties in the three regimes at (n) = 1. The algo- 
rithm used in Ref. [|, [LI] allow us to calculate real-time dynamical responses accurately, 
including the optical conductivity a(u > 0), since all the eigenvectors in the fermionic 
sector for a given spin and phonon configuration are obtained exactly. From the sum- 
rule, e g kinetic-energy, and the integral of a{ui > 0), the uj = Drude- weight D w can 
be obtained. Carrying out this calculation it was observed that Dy/ vanishes at A c i 
signaling a metal-insulator transition (MIT). Here the insulating phase is spin-FM and 
orbital-AF, while the metallic one is spin-FM and orbital-disordered [33]. The density 
of states (DOS) for A > A c i was also calculated in Ref. |TTJ and it presents a clear gap 
at the Fermi level. The qualitative shape of D\y vs A on 4 2 and 4 3 clusters was found 
to be the same as in ID and, thus, it is reasonable to assume that the MIT exists also 
in all dimensions of interest. 

Consider now the influence of hole doping on the (n) = 1.0 phase diagram, with 
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special emphasis on the stability of other densities as fi is varied. Fig. 6 shows (n) vs 
/j, in the intermediate- A regime. It is remarkable that two regions of unstable densities 
exist at low-T (similar conclusions were reached using the Maxwell's construction J53]). 
These instabilities signal the existence of PS in the Hkjt model. At low-density there 
is separation between an (i) empty e 9 -electron band with AF-ordered t2 9 -spins and a (ii) 
metallic spin-FM orbital uniformly ordered phase. This regime of PS is analogous to the 
spin-driven PS found at low-density in the 1-orbital problem [|K]]. On the other hand, 
the unstable region near (n) = 1.0 is not contained in the 1-orbital case. Here PS is 
between the phase (ii) mentioned above, and (iii) the insulating intermediate-A spin-FM 
and orbital- AF phase described in Fig. 5 ||36|| . In 2D the results were found to be very 
similar (Ref. (JTTJ] ) . The driving force for this novel regime of PS are the orbital degree 
of freedom, since the spins are uniformly ordered in both phases involved. Studying 
(n) vs /i, for A < A c i only PS at small densities is observed, while for A > A c2 the PS 
close to (n) = 1 is similar to the same phenomenon observed in the 1-orbital problem 
since it involves a spin-AF orbital-FM phase 0. 




Figure 6: (n) vs /i at A = 1.5, L = 10, and T = 1/40 (solid circles). The discontinuities 
near (n) = 1.0 and 0.0 show the existence of unstable densities. The solid line is 
obtained from the Maxwell's construction. The triangles are results also at A = 1.5 
and T = 1/40, but using 14 sites and only 2 x 10 4 MC sweeps to show the appearance 
of hysteresis loops as in a first-order transition. The inset shows the T-dependence of 
the results at L = 10. These results were taken from Ref. [11]. 

The phase diagram of the ID Hkjt model is given in Fig. 7. The two PS regimes 
are shown, together with the metallic spin-FM region. This phase is separated into two 
regions, one ferro-orbital ordered and the other orbitally disordered. The existence of 
these two regimes can be deduced from the behavior of the pseudospin correlations, the 
mean value of the pseudospin operators, and the probability of double occupancy of 
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2s 



the same site with different orbitals. The results are similar for several {tab} sets 
The simulations of Ref. |Tl| suggest that the qualitative shape of Fig. 7 should be valid 
also in D = 2 and 3. Recently we learned of experiments where the low-T coexistence 
of domains similar to those described in Ref. [11], i.e. one orbitally-ordered and the 



other FM-metallic, has been observed in Mn-perovskites |37L 



S-AF 

O-F 

insulator 
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O-AF 

insulator 
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Figure 7: Phase diagram of the H KJT model at low-T, J H = 8.0, J' = 0.05, and 
using set Ti for {t a b}. S — F and S — AF denote regimes with FM- and AF-spin orders, 
respectively. O — D, O — F, and O — AF represent states with disordered, uniform and 
staggered orbital orders, respectively. PS means phase separation. These results were 
taken from Ref. [11]. 



Optical Conductivity 

In Ref. []T| results for a(u) were also presented. Experimental studies for a variety 



of manganites such as Nd .7Sr 3 Mn0 3 , Lao.7Ca .3Mn0 3 , and Lao.7Sr .3Mn0 3 reported 
a broad peak at to ~ leV (for hole doping x > 0.2 and T > T C FM ) [10]. At room- 
T there is negligible weight near u = 0, but as T is reduced the leV-peak shifts to 
smaller energies, gradually transforming into a Drude response well below T c . The 
finite-u; peak can be identified even inside the FM phase. The coherent spectral weight 
is only a small fraction of the total. Other features at larger energies ~ 3eV involve 
transitions between the J#-split bands and the O-ions. In addition, Jung et al. EC 



interpreted the leV feature at room-T as composed of two peaks due to intra- and 
inter-atomic transitions in JT-distorted environments. 

In Fig. 8, ct{uj) for the H^jt model is shown at (n) = 0.7 and several tempera- 
tures near the unstable PS region of Fig. 7 (weight due to Jh split bands occurs at 
higher energy). Here the FM spin correlation length grows rapidly with the lattice 
size for T* < 0.05£, which can be considered as the "critical" temperature. Both at 
high- and intermediate-T a broad peak is observed at uj ~ 1, smoothly evolving to 
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lower energies as T decreases. The peak can be identified well-below T* as in exper- 
iments |40f| . Eventually at very low-T, <y(oj) is dominated by a Drude peak. The 
T-dependence shown in the figure is achieved at this A and (n) by a combination of a 
finite-d; phonon-induced broad feature that looses weight, and a Drude response that 
grows as T decreases (for smaller As, the two peaks can be distinguished even at low-T) . 
The similarity with experiments suggests that real manganites may have couplings close 
to an unstable region in parameter space. In the inset, Dyy vs T is shown. Note that at 
T ~ T*, Dy/ vanishes suggesting a MIT, probably due to magneto polaron localization. 
Results for the 1-orbital case are smoother, with no indications of a singularity. These 
features are in agreement with experiments, since the manganite "normal" state is an 
insulator. Work is currently in progress to analyze in more detail this MIT transition. 




Figure 8: (a) <j(uj) parametric with T, at A = 1.0, (n) = 0.7, and L = 20. The inset 
shows Dy\r vs T for both the Hkjt (soild circles) and the 1-orbital model of Ref. [8] 
(open squares) (the latter at (n) = 0.65). Dw is normalized to its maximum value at 
T = 0.01; (b) <j(oj) vs u parametric with (n) at A = 1.5, T = 1/10, and L = 16 (results 
for L — 10 are very similar). The inset shows the lower J^-split DOS at (n) = 0.93. 
Both in (a) and (b) a ^-function broadening e = 0.25 was used. These results were 
taken from Ref. [11]. 



A similar good agreement with experiments was observed working in the regime 
of the orbitally-induced PS but at a temperature above its critical value Tps (roughly 
~ 1/20, see inset Fig.6). Here the broad feature observed at high-T in Fig.8a moves 
to higher energies (Fig. 8b) since A has increased. At the temperature of the plot the 
system is an insulator at (n) = 1, but as hole carriers are added a second peak at lower 
energies develops, in addition to a weak Drude peak (which carries, e.g., just 1% of the 
total weight at (n) = 0.61). This feature at high-T is reminiscent of recent experimental 
results ||0| where a two-peak structure was observed at room-T and several densities. 
Similar results were obtained on 4 2 clusters. In Fig. 8b the peak at large-u; is caused 
by phononic effects since its position was found to grow rapidly with A. It corresponds 
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to intersite transitions between Mn 3+ JT-split states, 
compatible with a Mn 3+ -Mn 4+ transition |4l[ 



The lower energy structure is 
The inset of Fig.8b shows the DOS of 
the system. The two peaks above \i are responsible for the features found in o[uS). This 
interpretation is the same as given in Ref. [p7| at D = oo. 



EXPERIMENTAL CONSEQUENCES OF PHASE SEPARATION 

A large number of papers in the context of experimental studies of manganites have 
reported the presence of some sort of inhomogeneity in the system that is tempting 
to associate with the phase-separation tendencies observed in our studies. In the next 
paragraph we will mention some of these results such that the interested reader can have 
at hand at least part of the relevant references on the subject to decide by him/herself 
on this issue. 

Some of the experiments that have reported results that could be compatible with 
ours are the following (in a random order): 1. De Teresa et al. studied Lai_ x Ca x Mn03 
using small-angle neutron scattering, magnetic susceptibility and other techniques, at 
x=l/3. The analysis of their data showed the existence of "magnetic clusters" of size 



approximately 12 A above the ferromagnetic critical temperature; 2. Hennion et al. |43 
recently presented elastic neutron scattering results below T c at x = 0.05 and 0.08 also 
for Lai_ x Ca x Mn03. They interpreted their results as indicative of "magnetic droplets". 
The density of these droplets was found to be much smaller than the density of holes 
implying that each droplet contains several holes; 3. Lynn and collaborators [H have 



studied Lai_ x Ca x Mn03 at x=l/3 also using neutron scattering. Lattice anomalies and 
magnetic irreversibilities near T c were interpreted as evidence of two coexisting distinct 
phases; 4. Perring et al. ^TJ studying Lai.2Sri. 8 Mn 2 07 with neutron scattering reported 



the presence of long-lived antiferromagnetic clusters coexisting with ferromagnetic crit- 
ical fluctuations above T c ; 5. Allodi et al. [[45| using NMR applied to La!_ x Ca x Mn03 



at x ~ 0.1 observed coexisting resonances corresponding to FM and AF domains; 6. 



Cox et al. [46] studying Pr .7Cao.3Mn03 with x-ray and powder neutron scattering 



reported the presence of ferromagnetic clusters; 7. Bao et al. f|7| in their analysis 



of Sr 2 _ x La x Mn04 (2D material) found phase separation at small e 9 -densities; 8. Ya- 
mada et al. |^8| in their study of Lai_ x Sr x Mn03 with neutron scattering at x=0.1 and 



0.15 interpreted their results as corresponding to polaron ordering; 9. Roy et al. 
studying Lai_ x Ca x Mn03 near x = 0.50 found the coexistence of two carrier types: 
nearly localized carries in the charge-ordered state and a parasitic population of free 
carriers tunable by stoichiometry; 10. Booth et al. p0[ also studying Lai_ x Ca x Mn03 



found that the number of delocalized holes ndh in the ferromagnetic phase changes 
with the magnetization M as ln(ridh) oc —M. Since the magnetization saturates well 
below T c , this implies that there is a range of temperatures where sizable fractions of 
localized and delocalized holes coexist; 11. Jaime et al. |51| reported the possibility 



of polaronic distortions of the paramagnetic phase of La — Ca — Mn — O manganites 
persisting into the ferromagnetic phase, and analyzed the data with a two-fluid model; 
12. Zhou and Goodenough |52| studying the thermopower and resistivity of 16 0/ 18 
isotope-exchanged (Lai_ x Nd x )o.7Cao.3Mn03 found indications of a segregation of hole- 



rich clusters within a hole-poor matrix in the paramagnetic state; 13. Billinge et al. |53 
reported the coexistence of localized and delocalized carriers in a wide range of densi- 
ties and temperatures below T c ; 14. Ibarra et al. |54j] studied (Lao.5Ndo.5)2/3Ca 1 / 3 Mn0 3 



concluding that insulating charge-ordered and metallic ferromagnetic regions coexist at 



low temperatures; 15. Rhyne et al. p5| in their study of Lao.53Cao.47Mn03 with neutron 
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diffraction found two magnetic phases below the transition temperature, one ferromag- 
netic and the second antiferromagnetic, both persisting down to 10K; 16. Heffner 
et al. |)(J studying Lao.67Cao.33Mn03 using zero-field muon spin relaxation and resis- 
tivity techniques found indications of polarons on the spin and charge dynamics; 17. 
Jung et al. []57J have very recently studied the optical conductivities of Lay/gSrx/gMnOa 
concluding that there are indications in this system of phase separation. 

The common theme of these experiments, and possibly others that have escaped our 
attention, is that some sort of inhomogeneity appears in the analysis of the data. The 
phenomenon is apparently more clear at small hole densities x and low temperatures 
which is precisely the region where our results indicate a strong tendency to phase 
separate. But also above T c in the interesting densities around x=l/3 magnetic clusters 
have been reported. It is nothing but natural to imagine a "smooth" connection between 
x ~ 0.1 and low temperature with x ~ 1/3 and temperatures above T c . In this case 
the large droplets found by Hennion et al. could have evolved as to become the 
magnetic clusters of De Teresa et al. |4!| by a reduction of their size and number of 
holes inside. We strongly encourage experimental work that could contribute to the 
answer of this or other important issues around the common theme of possible charge 
seggregation in real manganites. 



SUMMARY 



Summarizing, a comprehensive computational study of models for manganites have 
found that the expected double-exchange induced strong tendencies to ferromagnetic 
correlations at low temperatures are in competition with a regime of "phase separation" . 
This regime was identified in all dimensions of interest, using one and two orbitals (the 
latter with Jahn- Teller phonons), and both with classical and quantum localized t 2g 
spins. It also appears in the presence of on-site Coulomb interactions. This robustness 
of our results suggests that phase separation may also be present in real manganites. In 
the previous section experimental literature that have reported some form of charge in- 
homogeneity in the context of the manganites has been briefly reviewed. It is concluded 
that theory and experiments seem to be in qualitative agreement and phase separation 
tendencies (which may correspond to the formation of magnetic droplets or even stripes 
once Coulomb interactions beyond the on-site term are included in the analysis) should 
be taken seriously. They may even be responsible for the phenomenon of Colossal Mag- 
netoresistance that motivated the current enormous interest in the study of manganites 
in the first place!. 
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